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Abstract 

In a paper entitled Beables for Quantum Field Theory, John Bell has shown that it was possible 
to build a realistic interpretation of any hamiltonian lattice quantum field theory involving Fermi 
fields. His model was constructed along the ideas he used to present the de Broglie-Bohm pilot 
wave theory. However, the beable (or element of reality) is now the fermion number density, which 
is not a particle density, as in the de Broglie-Bohm pilot wave theory. The model is stochastic but 
Bell thought that it would become deterministic in the continuum limit. We show that it is indeed 
the case, under an assumption about the physical state of the universe, which follows naturally 
from the Bell model. Moreover, the continuum model can be established directly. The assumption 
is that the universe is in a state obtained from the positronic sea (all positron states occupied) 
by creating a finite number of negative charges. The physical interpretation is the following: the 
negative charges are in motion in the positronic sea and their positions are the beables of the 
Bell model. The velocity laws we obtain for the motion of the negative charges are very similar 
to those given by Bohm and his co-workers for free relativistic fermions (first quantization). The 
Bell model is non-local (it is unavoidable); we show it explicitly in the simplest case. Under the 
previous assumption about the state of the universe, and for quantum field theories involving only 
Fermi fields, wave functions can be defined, and calculations can be performed as in non-relativistic 
quantum mechanics, since we stay in a sector of the Fock space with a fixed fermion number. 
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I. INTRODUCTION 

At the time Bell wrote his paper jy, Bohm had already shown that it was possible to 
build a realistic interpretation of any bosonic quantum field theory Q]. To achieve that goal, 
Bohm took the field as the beable (or element of reality), however he was not able to do 
the same for fermions. The aim of Bell was then to show that it was also possible to build 
a realistic interpretation of any fermionic quantum field theory, along the pilot-wave ideas 
given by de Broglie and later by Bohm. Bell managed doing so but he took a really different 
beable: the fermion number density. It is quite different from the non-relativistic pilot-wave 
theory, whose beables are the positions of the particles. The model is also formulated on 
a spatial lattice (space is discrete but time remains continuous). His model is stochastic, 
but he suspected that the theory would become deterministic in the continuum limit. The 
crucial point is that the fermion number density is not a particle density. It is in fact related 
to the charge density. 

To deal with fermions on a spatial lattice, we have used the Banks-Susskind-Kogut theory 
of staggered fermions It is the best way to eliminate the fermion doubling problem. 
Taking the continuum limit of the Bell model for a staggered lattice, we show that the theory 
becomes indeed deterministic in that limit. The velocity laws we obtain are very similar to 
those given by Bohm and co-workers for free relativistic fermions (first quantization). We 
have worked with a one-dimensional lattice to simplify the expressions but we think that it 
could be extended to the three-dimensional case without difficulty. This remark is justified 
since we show that the continuum model can also be established directly. 

Finally, we will study in more details the fact that any purely fermionic quantum field 
theory can be brought to a particular sector of the Fock space, with a fixed fermion number, 
allowing us to define wave functions and to perform calculations in the same way as they 
are done in non-relativistic quantum mechanics. We will also study entanglement and non- 
locality in the simplest case (two quanta of the fermion number). 
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II. THE BELL MODEL 



A. Ontology 

Three-space continuum is replaced by a finite lattice, whose sites are labelled by an index 

/ = 1,2,...,L. 
The fermion number density is the operator 

F(i) = ^{i)m = 4(i)Mi) + ■■■ + 4(i)Mi) ■ 

Since 

[F(k),F(l)]=0 Vfc,Z E {1,2,...,L}, (1) 
it is possible to define eigenstates of the fermion number density: 

F(l)\n,q} = f(l)\n,q) , 

where q are eigenvalues of observables Q such that {F(l), F(2), . . . , F(L), Q} is a 
complete set of observables, and n is a fermion number density configuration (n = 
{/(l), /(2), . . . , f(L)}). Eigenvalues f(l) belong to {0, 1, 2, 3, 4}. 

Thus we can imagine that the universe is in a definite fermion number density configu- 
ration n(t) at each time t, and that a measurement of the fermion number density at time t 
would simply reveal the configuration n(t). In Bell's words, the fermion number density is 
given the beable status. 

The second element, in the description of the universe, is the pilot-state ^(t)). Hence 
the universe, at time t, is completely described by the couple (\^(t)),n(t)). 



B. Equations of motion 

For the pilot-state, the Schrodinger equation is retained: 




= H\*(t)) . 



An equation of motion for n(t) must be added (it is called the velocity-law in pilot- wave 
theories). Call P m (t) the probability for the universe to be in configuration m at time t. 



3 



Then the velocity-law must be such that the relation 

P m (t) = ^|(m,g|vl/(t))| 2 (2) 

q 

holds for any time t, in order to reproduce the predictions of orthodox quantum field theory. 
Since the configuration space is discrete, it is impossible to find a deterministic velocity- 
law. Instead jump-rates have to be defined. Call T nm (t) the jump-rate for the transition 
m — > n at time t, for m 7^ n (T mm (t) = Vm). In other words, T nm (t) is the probability 
density (probability by unit of time) for the universe to jump in configuration n, knowing 
that the universe is in configuration m at time t. The probability for the universe to stay in 
configuration m at time t + dt is obtained by the normalization condition and is equal to 

l-J2 T nm(t)dt , (3) 

n 

for dt small enough. Eq. (J2J) is assumed to be true for an initial time to; then the constraint 
on the stochastic velocity-law becomes 

^=E|ik 9 i*w>i 2 - w 

q 

Let's calculate the first member: 

P m (t + dt) = J2 T mn {t)P n {t)dt + (1 - T nm (t))P m (t)dt , (5) 

n n 

from which follows 

= Jim Fm(t + ^~ Pm(t) = £(7U*)P n (t) - T nm (t)P m (t)) . (6) 

n 

Now let's calculate the second member of equation Q. With the help of the Schrodinger 
equation, we have 

d\(m,q^(t))\ 2 ^ mjg |^ f )^( £ )| ii7 | m) ^ + ( m!g | -iH\V(t))(*(t)\m,q) (7) 

=29te[(m,q\ - iH\V(t))(V(t)\m,q)] (8) 
=2^2^Rt[{^(t)\m,q)(m,q\ - iH\n,p)(n,p\V(t))] . (9) 

n,p 

Hence the constraint on the velocity law is 
J2( T mn(t)Pn(t) ~ T nm (t)P m (t)) = 2 ^ <Ke [(*(*) |m, q) (m, q\ - iH\n, p) (n, p\ *(*)>] , (10) 



or 



T mn (t)P n (t) -T nm (t)P m (t) =2j2M(V(t)\m,q)(m,q\-iH\n, P )(n,p\V(t))] . (11) 



If one takes for the following definition for the jump-rates: 

Tmn(t) = J ~J^ if Jmn > , (12) 

T mn (t) = otherwise , (13) 

where 

J mn (t) = 2j2Xt[(y(t)\m,q)(m,q\-iH\n,p)(n,p\y(t))\ , (14) 



p,q 



then Eq. (jlj) is satisfied. 
C. Comments 

We would like to make some remarks about the fermion number density. First, eigenstates 
of the fermion number density are also eigenstates of the charge density, which is 

-e : ^(l)m : , 

at least if we consider only electrons and positrons. The fermion number 

s=2 „ 

F =J2 l C l(P)C s (p) + d s (p)dl(P)} , 

s=l ^ 

is not the particle number, which is 

s=2 . 

N = J2 d 3 p [ct(p)c s (p) + d\(p)ds(p)\ ■ 

s=l ^ 

In fact the fermion number density does not commute with the particle number; it is possible 
to find well-behaved functions f(x) such that 

d 3 xf(x)^(x)ifj(x),N] . 

The proof is given in appendix El 

Thus the charge density does not commute with the particle number either. A measure- 
ment of the charge contained in any finite region including the coordinate Xq, with value— e, 



is never an electron. In fact, as we will see, it is a superposition containing one electron, 
two electrons and one positron, three electrons and two positrons, and so on. On one hand, 
it is disturbing, since the tracks observed in bubble chambers are said to represent electrons 
or positrons. On the other hand, when the electromagnetic field is taken into account, it is 
quite natural, since measurement of localized properties involve high energy radiation, and 
thus that can lead to pair creations. But these are are just few remarks to draw attention 
to the interpretation problems one has to cope with. 

III. THE DIRAC THEORY IN A 1 + 1 SPACE-TIME 
A. The Dirac equation 

1. Solutions of the Dirac equation 
The hamiltonian is 

H = ap + [3m. (15) 
Since W = H and H 2 = p 2 + m 2 , we obtain the following relations: 

a f = a ft = {3 {a, (3} = a 2 = (3 2 = 1 . (16) 

The smallest dimension for a representation of that algebra is two. For example: 

(o 1 




a = 



1 



(17) 



Thus we have a spinor of dimension 2, if)(t, x), which is solution of 



i W,*) = _ ia Wt,z) 

at ox 

This equation can be rewritten in a covariant form, by introducing the 7 matrices, defined 
by 

7° = (3 = 70 7 1 = a = -71 . (19) 

That implies 

(^<9 M -m)i(j(t,x) = . (20) 
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One can verify that the 7 matrices satisfy the following relations: 

{7", 7"} = 20"" 7^t = yy 7 o _ (21) 

Each component of the spinor is a solution of the Klein-Gordon equation: 

(n + m 2 )ip(t,x) = . (22) 
Now we search for free solutions; the most general form is thus 

u(p)e- lE ^e ipx v(p)e tEpt e~ tpx , (23) 

with p e and E p = ^Jp 2 + m 2 . We find that 

il) p + (t,x) = u(p)e~ iEpt e ipx = I 1 e-"»e** (24) 

\m+E p J 

ip p -(t,x) = u{p)e iEpt e- ipx = ^™+E^j e iE p t e -i P x ( 25 j 

are solutions of Dirac equation, respectively of positive and negative energy. To obtain an 
interpretation of the theory, we need a conserved current, whose temporal component is 
positive: 

d^f = with f > . (26) 

The current 

f = ^rv , (27) 

with Tp = "0^7°, is suitable. Spinors normalization still remains to be discussed. Assume 
that the universe is a box of volume V, in a inertial frame S where the momentum of the 
free particle is p. Then the quantity 



/ dxf{t,x) 
Jv 



must be equal to 1 in every inertial frame. That means that the spinors u(p) and v(p) must 
be normalized to ^r-^f, where V is the volume of the universe in an inertial frame where 
the particle is at rest. 
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2. Physical interpretation 

Dirac obtained a conserved current, whose time- component is positive, a task that was 
impossible with the Klein-Gordon equation. But the negative energy states are still there. 
Once interactions are taken into account, that would lead to the instability of the hydrogen 
atom, for example. To avoid this, Dirac assumed that all negative energy states were occu- 
pied. Hence a positive energy electron cannot transit to a negative energy state, due to Pauli 
exclusion principle. That state of lowest energy is called the Dirac sea. It is impossible to 
distinguish it from a state where no electrons are present. The absence of a negative energy 
state of momentum p ( a hole in the Dirac sea) would be seen as a particle of positive energy 
a/p 2 + m 2 , momentum — p and charge e. That led to the prediction of anti-particles known 
as positrons. 



B. The Dirac quantum field theory 

The first step, in the construction of the corresponding quantum field theory, is to obtain 
a classical relativistic action, from which we can obtain the Dirac equation, by using the 
least action principle. The following action is suitable: 

^-/^t,.)(*^-» W «,.). (28, 
It is not hermitian but it can be rewritten as 



s = s h + J «v M > 



where Sh is hermitian, and the last term can be dropped. We get the momenta conjugate 
to the fields: 

7r a (t,x)= ff*; 30 s) <(t,x)= ffff; g) =o . (29) 

d(il> a (t,x)) d(^(t,x)) 
The next step is quantization, according to the canonical equal-time anti-commutation re- 
lations: classical fields become quantum fields, obeying the relation 

{ip a (t,x),n b (t,y)} = i8(x - y)S ab {ip a (t,x),ifo(t,y)} = {ir a (t, x), n b (t, y)} = . 

Those relations can be rewritten as 

{4> a (t,x),4>l(t,y)} = 5{x-y)5 ab {ip a {t,x)^ b {t,y)} = . (30) 
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Since the quantum field if)(t, x) satisfies the Dirac equation, it is a superposition of free 
solutions with operators as coefficients: 

if>(t,x) = -tL / dp[c{p)u{p)e- iEpt e ilpx + C{-p)v{p)e lEpt e- ipx ] , (31) 



4>\t, x) = -= / dp[c^{p)u T (p)e iEpt e- ipx + C{-p)v T {p)e- iEpt e ipx ] . (32) 

V27T J 

c, C, and are operators satisfying unknown anti-commutation relations, that must be 
chosen in order to regain the equal-time anti-commutation relations (Eq. (|3UJl ). Spinors are 
normalized to 

uHp)u(p) = — vHp)v(p) = — 

m m 

With the help of Eq. (]3*T]) and Eq. (|32|1 . we can work out the equal-time anti-commutation 

relations. We have 

{Mt,x),Mt,v)} = ^f dpdg[{c(p),c(g)}e-^ e - iB » t e* !e e < «« a (p) U6 (g) 

+{a-p),a-q)}e lEpt e lE ^e- ipx e-^v a (p)v b (q) 
+{c(p),a-q)}e- iEpt e lE ^e ipx e-^u a (p)v b (q) 
+{({-p),c(q)}e iE v t e- iE * t e- ipx e iqy v a {p)u b (q)] . 

If we take 

{c(p),c(q)} = { c (p) 1 ((q)} = {C(p), ({q)} = Vp,?, (33) 

then we obtain {^(t, x), ^(i, ?/)} = 0. The relation {-^ a (t, x), ^(t, y)} = 5(x — y)5 a 6 remains 
to be considered. With the help of Eqs (J3T|) and ([32]) . we have 

{lM*.*),W(*>v)} = ^| ^g[{c(p),ct(g)}e-^e^ t e^e-^ a (p)n fe (g) 

+{C(-p),C t (-g)}e^ i e-^ 4 e-^e^t; a b)^(g) 
+{c(p),C t (-g)}e" i£:p *e-^*e^e^M a (p)^(g) 
+{C(-P), c+C^le^e^e^e-^Ua^^Cg)] . 

Taking 

777 777 

{c(p), c t( g )} = — tf(p- g ) {CCpU'Gz)} = Wp),C t (?)} = o v P , g , 



the relations (|5U|) are regained. We can also choose 

x) = -L J dp^-[c{p)u{p)e^ Et e^ + C(- p )v{p)e iEt e^} , (34) 

with the following anti-commutation relations: 

{c(p), c(q)} = {c(p), ((<?)} = {C(p), C(?)} = 

{c(p), c t(g)} = 5(p - g) (C(p), C f (g)} = *(p - g) {c(p), C f (<?)} = V p, q . 

That is the choice we adopt. Now the observables can be expressed in the momentum space. 
For the hamiltonian, we have 

H = J dxtp\x)\— iaV + m/3]ip(x) = J dp\/p 2 + m 2 [c^ (p)c(p) — C^( — p)C(~p)] • 

The momentum is 

P = y cto/^(x)[— zVj^x) = y dpp[c^(p)c(p) — C^( — p)C( — p)] • 

And the fermion number is 

F= dxip^ (x)ip(x) = / dp[c\p)c{p) + C^( — p)((— p)] ■ 



We can define a vacuum as a state annihilated by any operator c(p) or C(p)i we can that 
state |0i): 

c(p)|0i) = C(p)|0i) = Vp. 



Now it is clear that (p) creates an electron of energy y/ 'p 2 + m 2 and momentum p, whereas 



C^(p) creates an electron of energy — a/p 2 + m 2 and momentum p, and that the Dirac sea is 

the state \DS) = TTC^(p)|0i}- Usually, everything is rewritten by introducing positrons, by 
v 

making the substitutions 



where d^(p) is the operator that creates a positron of momentum p and energy a/p 2 + m 2 . 
And we have to define another vacuum, | O2) : 

c(p)|0 2 ) = d(p)|0 2 ) = Vp. 

We have the relations 

io 2 >=nc t (p)ioi> io 1 >=n^(p)io2> • 
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IV. FERMIONS ON A LATTICE 



One-space continuum is replaced by a lattice of spacing 5, having L = 2N sites. The 
momentum space is also a lattice, having 2N sites and a spacing j^. To get the lattice 
action, one makes the following substitutions 

dx/>(t,x) il>{t,j + l)-i/>(t,j-l) 



dx 25 
in the continuum action 

S[tp,^] = J dxdtip(t,x)[i'fd li -m]il}(t,x) . 

Doing so, we obtain 

We can eliminate the factor 5, by making the substitution 
so that the lattice Dirac action is 

The least action principle gives the lattice Dirac equation (we use periodic conditions on the 
boundaries = ipx)' 

whose free solutions are 

^(*,j) = «(Piat)e _i£,ot(p) V w '*= [ 1 |e- ,Ew(p) 'e iB ' { (35) 

Plat 



^m+Eiatip), 

Plat 



i/^(t,j)=u( P i at )e iEl ' tt e- ipj5 = \ ™ +E ^W J e iE lat (p)t e -i Pj s ) (36) 

where p/ at = sm ^ an d Ei at {p) = \Jpi at + "^ 2 - Now we turn to the quantization. The 
anti-commutation relations become 

{Mt,j),^ b (t,k)} = 5 a b 6i , 
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and all other anti-commutators vanishing. Since it satisfies the lattice Dirac equation, the 
quantum field ip(t,j) is a superposition of free solutions, with operators as coefficients: 

V 

Again, u(p) and the anti-commutation relations satisfied by the operators c, d, and <v , 
are determined by the canonical equal-time anti-commutation relations. But is not difficult 
to see that the hamiltonian is 

H = E y/^M +m 2 [c t (p)c(p) _ d(p)rft(p)] 
v 

Thus there are four states of energy m. Generally, every eigenstate containing n particles 
is degenerate, with degeneracy 2™. The same problem occur in the propagator; it has four 
poles and thus propagates twice more particles. It is called the fermion doubling problem 
and there are many theories to deal with it (the Wilson theory, the Banks-Susskind-Kogut 
theory of staggered fermions, to mention the main ones). In the continuum limit (lattice 
spacing going to zero and finite momentum), we have 

sin(p<5) 

so that the fermion doubling problem disappears. But that is not a reason to ignore it in 
our calculations. In the one- dimensional case, it is easier to use the Banks-Susskind-Kogut 
theory of staggered fermions to overcome it. 

A. The Banks-Susskind-Kogut theory of staggered fermions 

The idea is to start from the previous theory, with a lattice spacing equal to 25, and 
to say that there are two superposed lattices, one where upper components of ip(t,j) live, 
and another one where lower components live. By moving the lower lattice to the right, 
with a translation of magnitude 5, we obtain a theory of a complex field (f)(1), over a lattice 
containing twice more sites. 

Since the part of the article which we are interested in is quite small, we will just quote 

it 

Consider a spatial lattice (continuous time) with a lattice spacing a. Label the 
lattice sites with an integer n. There will be a one-component fermion field (f>(ri) 
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at each site n. 4>{n) satisfies the anti-commutation relation 

{tf(n), 0(m)} = 8 nm {0(n), 0(m)} = (37) 

</>(n) is related to a properly normalized continuum field x having canonical 
anti-commutation relations by 

(f>(n) = y/ax(x) (38) 

Consider the hamiltonian 

H = Ya X> f ( n M n + 1) - f (^ + 1)0(»)] (39) 

n 

We claim that with a proper identification of a two-component fermion fields 
Eqs ()37)) -([39 )) generate the massless Dirac equation in the continuum limit. First 
compute 

i[H, <f>(n)] = 0(n) = ^ ^ ^ (40) 

Note that the time dependence of 4>{n) at even (odd) sites is determined by the 
spatial difference of 0(n± 1) at odd (even) sites. So, to ensure finite time depen- 
dence in 4>(n) at even (odd) sites, we must require that the spatial dependence 
in 0(n) at odd (even) sites be smooth. Thus we defined a two-component field 
ij)(n) as follows: 

■0e 




ip e (n) = 0(n), n even (41) 
V'o(^) = n odd. 

Then the components of ip{ri) satisfy the equations 

= ~Ax' = "Ax" (42) 



where A indicates the discrete difference in Eq. fiT)|) . Note that Eq. (|4^|) 
becomes the massless Dirac equation in the continuum limit, 

d (o l\ d 

*H 10 (43) 



13 




in a standard basis where 

7o 



If we consider the case of massive fermions, we just add a term m ^^( — l) n 0^(n)0(n) to the 

n 

massless hamiltonian: 

H = ~Y 5 W0(n + 1) - (j)\n + l)4>{n)\ + ^ m(-l) n tf(n)<j>(n) . 

n n 

There is a minus sign in the kinetic term, compared to the expression in the Banks- Susskind- 
Kogut article, but that is just a matter of conventions. Again we take the number of lattice 
sites to be 2N. Thus the momentum space is a lattice containing 2N sites, with spacing 

Slllf T)0 ) 

What about the fermion doubling problem? The momentum lattice is still — , but the 

o 

momentum space is divided into two sub-spaces: 

• the sub-space V±, where electrons live. 

• the sub-space V2, where positrons live. 

a , • • • •! 1 f , ^ , ^ x 11 sin(rx)) sin(g5) . 

And it is impossible to find p, q e V\ (resp. Vi) such that — ^ — = — . It is quite 



straightforward to show if we think in terms of degrees of freedom. Thus there are iV sites 
in V\ and iV sites in 7^2- 



1. Eigenstates of the fermion number density 

The fermion number density is the following operator: 

<t>\im) ■ 

The fermion number density is an operator with positive eigenvalues. Summing over all 
sites, we obtain the fermion number: 

1 peVj peP 2 

The eigenvalues of the fermion number range from to 2N. The state with the lowest 
fermion number is thus the positronic sea (all positron states occupied): 

ioi> = n ^(^)i°2> • 

peP2 
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Since the eigenvalues of the fermion number density are positive, there can be only one 
eigenstate of the fermion number density with fermion number equal to 0; it is thus the 
positronic sea: 

*t(ZM0|Oi> = o. 

From the anti-commutation relation 

we find that <jy(k) creates a quantum of the fermion number at site k, whereas 4>(k) destroys 
it. The eigenvalues of $ (I) (f)(1) belong to {0, 1}. The positronic sea is annihilated by any 
annihilator <fi(k) (it could be called the vacuum of the fermion number density). Eigenstates 
of the fermion number density are thus generated by applying creators at different sites. 

V. THE CONTINUUM LIMIT 

The fermion number commutes with the hamiltonian 

I 

and any physical state is an eigenstate of the fermion number, so there is an integer uj such 
that 

X>t(0tf(0l*(*)> = "!*(*)>■ 
i 

for any time t (with u G {0,1,..., 2N}). 
A. A localized quantum 

Let's first consider the case of a localized quantum 

i 

Since we are restricted to the sector of the Fock space with the fermion number equal to 1, 
we use a more convenient notation for the eigenstates of (f>\l)(f){l) 

|*>=0 t (fc)|O 1 > <p(l)<Kl)\k) = 6 lk \k) 
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The pilot-state can be decomposed in the basis {|&)}: 

k=N-l 

\*(t))= Yl ■ 

k=-N 

Assume that the beable at the initial time is n/(£o) = 8i k , corresponding to the state \k) = 
(p^(k)\0i). How does it evolve with time? The hamiltonian matrix element which appears 
in the transition current (Eq. iTlj) is 

(l\ - zH\k) = (Oil - i^H^ik)^) = ^ - ^ . 

Thus they are only two transitions which are allowed: 

J ik+ i )k = ^ito[V*(t,k + l)V(t,k)] 

J {k -i )k = -^Wt[**(t,k-l)*(t,k)] 

The quantum can only jump to the first neighbor site, either to the right or the left. Using 
the relation 

k - 1) = k + 1) - 2<5Vtf *(t, k) 
and keeping only leading terms, we get: 

J (fc+ i )fc = i«e[**(t,fc+l)*(t,A;)] 

J {k -i)k = -^9te[tf*(t, k + (t, fc)] . 

These two currents have opposite signs; then there is only one transition. The theory is thus 
deterministic. Let's consider the case 9ie[&*(t, k + l)\I/(£, k)] > 0. Then the quantum moves 
towards the right and its velocity is 

= ST (k+1)k dt = 5J (k+1)k = Vte[V*(t,k + l)V(t,k)] 
V ~ dt ~ \V(t,k)\ 2 ~ \*(t,k)\ 2 1 } 

Now what about the continuum limit? Remember that even sites correspond to upper com- 
ponents of the spinor, whereas odd sites correspond to lower components of the spinor. In the 
continuum limit, the couple (0t(n), ^(n+l)), for n even, tends to the spinor (ip\(x), ip\(x)). 
There are two sites that converge to a given x. If k was even, we still have to consider its 
companion k + 1. Then let's assume that the beable is <W+i) at time t . One can find that 
the currents are 

J(k+2)k+i = -^mt[y*(t, k)*(t, k + 1)] 

J(k)k+i = ^Ke[**(t,fc)tt(t,A; + l)] . 
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We still consider the case $He[^*(t, k + 1)\P(£, &)] > 0. Then only the current J(k+2)k+i 
survives and we find that 

J( k+ i)k = J( k+ 2)k+i = ^9te[tt*(t, k + (*, fc)] . 

The continuum current is just the continuum version of the sum of the two lattice currents 
(the one for k plus the other one for k + 1): 

J(M) =2«e[(*(*)|^(t,a;)|0 1 )(0 1 |Vi(*,a;)|*(t)>] 
= (M/(t)|^ t (t,x)|0 1 )a(0 1 |^(t,x)|M/(t)) 

B. Two localized quanta 

Now let's consider the case of two localized quanta: 

i 

We use the same notations: an eigenstate of ft (l)4>(l) with eigenvalue n\ = 5i kl + Sik 2 is 
denoted 

|A; 1 ,A; 2 ) = ^(A; 1 )0t( fc2 )|o 1 ) 

with ki < k 2 and k\,k 2 G {— AT, — AT + 1, . . . , N — 1}, corresponding to the state 
(ft (ki)(ft (k 2 )\Qi) . The most general beable at the initial time is thus ni(t ) = 8i kl + 5i k2 . 
The hamiltonian matrix element which appears in the transition current is 

(h, i 2 \ - iH\h, k 2 ) = - XfoiMWtW^O — + 1) 2 ~/ (; ~ - ^(fci)^(fc2)10i) • 

I 

Using the canonical anti-commutation relations and the fact that the positronic sea is anni- 
hilated by any <f>(n), and taking into account the constraints k\ < k 2 and l\ < l 2 , one finds 
that 

/III TT\1 1 \ _ ^2(fc2 + l)^lfel + ^l(fel+l)^2fc 2 _ ^ 2 (fe2-l)^lfcl _ ^l(fel-l)^2fc 2 

f 2| _ % tl \Ki, K 2 ) — — , 

only if k 2 ^ k\ + 1, otherwise 

^ 2 (fe 2 +i)^ifci ~~ ^i(fei-i)^ 2 fe 2 



(Zi, l 2 \ - iH\ki, k 2 ) = 
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The pilot-state can be decomposed in the basis '(ni)<ft (n 2 )|0i)}: 

W)> = 5^ ^(ni,n 2 ,t)\ ni ,n 2 ) . 
The transition currents are thus 

•/(ki.taMfci-i.fca) = -^9le[tf*(fci - l,k 2 )^(k 1 ,k 2 )} 

J( kl M)^k 1+ iM) = ^e[V*(ki + 1, hWh, k 2 )\ 

J (k u k 2 )^(kiM-i) = -^9te[tf*(fci, k 2 - l)*(fci, k 2 )\ 

^(fci,*a)-»(fci,*a+l) = ^9^[**(fcl> ^2 + l)*(fci, Jfe)] 
whether fc 2 = £4 + 1 or not, since ni) = 0. Keeping only leading terms, we get 

J (kiM)->(ki-iM) = -^9te[**(fci + l,fc 2 )tf(fci, fa)] = -J(fe 1 ,fe 2 )-,(fe 1 +i,fe 2 ) 

J^M^^te-i) = -^9*e[#*(£;i, £; 2 + l)*(fci, fca)] = -^(fci,fc 2 )-(fci,fc 2 +i) 

Thus there are only two transitions allowed, since the four transitions can be arranged by 
pairs whose currents have opposite signs. If k\ and k 2 are even and correspond to x\ and x 2 
in the continuum limit, the two transition currents are the two discrete components of the 
current 

\^SHe[<*(*)|V'I(a;i)^J(ar 2 )|0 1 ><0 1 |^ 1 (x 2 )^ 1 ( 3 : 1 )|*(t)>] y 
Again, it is thus deterministic, since the cases odd-odd and even-odd are similar. There are 
four couples that correspond to the same (xi, x 2 ) ((fci, k 2 ), that we have already considered, 
(ki + 1, k 2 ), ((hi, k 2 + 1), and ((ki + 1, k 2 + 1)). Summing their currents, we get 

/^;<*(0IV' t ^i)^(^)|0i>a(0 1 ^ a (x 2 )^(x 1 )|*(t)>\ 



J ee {Xi, x 2 , t) 



J(x 1 ,x 2 ,t) 



5](*(0IVl(a:i)^ t (s2)|0i>a(0i|^(x 2 )^(xi)|*(t)> 
V a / 



There is also another way to prove that the model is deterministic in the continuum limit. We 
take the continuum limit and then make a rotation in the configuration space of dimension 2, 
in order to align the axis X\ along a preferred direction, for example that of J ee (xi, x 2 ,t), if 
we consider that case. Then it can be shown that there are two transitions allowed, namely 
the two transitions along J ee {x 1 , x 2 ,t). The two currents have equal magnitude but opposite 
signs, so only one transition remains. Then the particle moves in a deterministic way and 
the scheme can be applied again for time t + dt, and so on. 
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C. Generalization 



Suppose we have u quanta 

Y,<t>\iw)m))="m)) ■ 
i 

and an initial beable n;(t ) = Sikj. + ^ik 2 + • • • + bik^- In the continuum limit, if Xj is the 
coordinate corresponding to kj, we expect to get the following current for the j-th coordinate: 

= X X • • • X **i...»i...»u.(*> ^i' • • • ' I «) Q «j»*»i...«..^( t ) xi, . . . , x w ) , 

where 

) = (0 1 |^(^)...^ Sl (x 1 )|^(t)) . 

It is just a refinement of the two-quanta case. In fact, there will be 2uj transitions allowed 
(to first neighbor sites). Those transitions can be arranged by pairs. In each pair, the two 
currents have equal magnitude but opposite signs. Thus there remains only u transitions 
(see Eq f)14jl ). Those uj currents are just the discrete components of a continuum current 
in a configuration space of dimension u. The deterministic character can also be proved by 
making a rotation in the configuration space of dimension u, as it was explained for the two 
quanta case. 

D. The continuum limit right from the start 

We have a physical state which evolves according to the Schrodinger equation and 

we know that there is an integer uj such that 

J d?x^{x)4){x)\^{t)) =w|tt(t)) • 

Thus can be decomposed along eigenstates of the fermion number density with fermion 

number equal to uj; those eigenstates are 

{ipl 1 (xi)...ipl w (xZ)\0 1 ) xu...,x w e R 3 ai,...,a w G {1,2,3,4}} . 

Thus 

j si=4 Suj=4 „ 

= ^X--X / ^•■■^^ S i... s j^^i,--- ) ^)^i 1 (^i)---^l(^)i0i) , 

si=l s w =l 
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where the wave function is antisymmetric. The universe at time t is described by a point 
in a configuration space of dimension 3u and by a pilot-state \^(t)). In the standard 
interpretation, the probability density to observe the universe in a configuration (x±, . . . , x w ) 
is 

si=4 s w =4 

Ptfr, . . . ,2 U ) = ■ ■ ■ E I WWl^l) • . .^l(^)|0i)| 2 , 

Sl = l So;=l 

and we have 

J dfxi . . . d 3 x U} p t (x 1 , ■ ■ ■ , x u ) = 1 . 
The time derivative gives a probability density current in the configuration space: 

7 n S\ 4 St^J 4 n 

— \ d 3 Xi . . . d 3 x u) p t (xi, . . . , Xy) = ^ . . . ^ / rf 3 fi . . . rf 3 f w 

Sl=l s w =l 

where we have switched to the Heisenberg picture. It can be simplified, knowing that 

.dip(t, x) 



dt 



—ia ■ Vip(t, x) + mf3ip(t, x) 



Terms containing (3 cannot contribute. Thus we obtain the following current for the j-th 
coordinate: 

s=4 si=4 Su=4 

^ 1 ,...,f w ,o = EE-E w°) i^i *) • • • ^ ^ *) • • • ^ *) i 01 ) 

S = l S\=l Suj = l 

oi SjS {0i\ip Su (xZ,t) . ..1p s (Xj,t) . . .ip Sl (xi,t)\^(t )) . 

Define A = (<x l5 . . . , a w ), Va G i? . If is the position of the universe in the configuration 
space, at time t, its velocity is thus 



J(X,t) 



X=X(t) 



Pt(X) 

VI. QUANTUM NON-LOCALITY 

The Bell model is non-local, but this is a necessary property of any realistic interpretation 
of quantum field theory, following the EPR paradox, Bell's inequality and experiments. We 
just want to show it explicitly. Consider the case of two negative charges moving in the 
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positronic sea. Then there is an interaction among these two charges, by the Pauli Principle; 
the wave function ^f SlS2 (t, x%, X2) is antisymmetric. The least entangled state, satisfying the 
antisymmetry requirement, is 

^ SlS2 (t,xi,x 2 ) = Xsi(t,Xi)® S2 (t,x 2 ) - ® Sl (t,xi)x S2 (t,x 2 ) . (45) 

The currents are 

J 1 (t,x 1 ,x 2 ) = ^ i S2 (t , xi , x 2 ) ^ 2s2 (t,xx, x 2 ) + % S2 (t,x 1 ,x 2 )^ 1S2 (t,x 1 ,x 2 )] 
J 2 (t, x u x 2 ) = (t, Xi, x 2 )^ si2 (t, Xt,X 2 ) + ty* si2 {t, x 1 ,x 2 )^ Sl i(t, x 1 ,x 2 )] 

si 

Substituting ^/ SlS2 (t, x±, x 2 ) by the right-hand part of Eq. (J45|) . we get, for J\\ 

Ji(t, x x , x 2 ) =(x*iX2 + X* 2 Xi)(t, ^i)(|$i| 2 + |$ 2 | 2 )(t, x 2 ) 
+ ($t$ 2 + ^$ 1 )(t,x 1 )(| X i| 2 + |x2| 2 )(t,x 2 ) 
~(xl®2 + X* 2 ®i){t, Xl){®lXl + $* 2 X2)(t, x 2 ) 
-(®\X2 + ® 2 Xl)(t, X\){XV*\ + X* 2 ^)(t, X 2 ) 

This general form is inconsistent with the existence of two real currents Jf and J® such 
that 

Jl(t,Xi,X 2 ) = Ji(t,Xi)Ji (t,x 2 ) , 

thus the model is clearly non-local: the velocity of one of the charges, at time t, depends 
on its position, on the pilot-state, as well as on the position of the other negative charge, at 
the same time. 



VII. QUANTUM FIELD THEORY IN A FIXED SECTOR OF THE FOCK SPACE 

It is interesting to note that quantum field theory calculations can be done in the same 
way as they are made in non-relativistic quantum mechanics, at least if we only consider 
fermions. Let's consider the following model, fermions interacting through a quartic term: 

H = J d 3 x (ip\x)[-ia ■ V + m/5]^(x) + g(ip\x)Pi/;(x)) 2 ) . 

Assume that there are two negative charges in the positronic sea: 
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with 

Then the pilot-state can be decomposed along eigenstates of the fermion number density, 
with fermion number equal to two: 

si=4s2=4 „ 

= y^^2 d 3 x l( l 3 x 2 ^ slS2 (t, x 1 ,x 2 )i;l l (x 1 )iPl 2 (x 2 )\0 1 ) . 

• si=l S 2 = l 

Substituting the pilot-state by the right-hand part of the previous equation in the 
Schrodinger equation, and projecting onto a state (xi)ipl 2 (x 2 )\0i) , we get 

i^* a i S2 (*,zi,f 2 ) =(/9^(x 1 ,f 2 )) SlS2 - (/3*(£ 2 ,£i)) S2Sl 

-i(a ■ V Sl ^ t (xi,x 2 )) slS2 + i(a ■ V a?2 ^t(f 2 , 

fi, f 2 ) / 9 T ) SlS2 <5(f 1 - f 2 ) - £ 2 , fi)/3 T ) S2Sl 5(fi - f 2 ) , 

where we have dropped an infinite constant. 

VIII. CONCLUSION 

We have obtained the continuum limit of the Bell model, for fermions living in a one- 
dimensional space, using a staggered lattice and we have also shown that we could build the 
continuum Bell model directly. Physically, it is a theory of negative charges moving in a 
positronic sea. There is an underlying assumption about the state of the universe, namely 
that it is an eigenstate of the fermion number (which is always true), with a finite eigenvalue. 
That follows naturally from the Bell model itself. Can one build a similar interpretation for 
the Klein-Gordon quantum field theory? It seems that the answer is no, for it is impossible 
to define a state annihilated by a charge creator in the Klein-Gordon quantum field theory. 
Another point worth mentioning is that the construction of the Bell model has nothing to 
do with the equation of motion being linear. We could use a Van der Waerden field and 
obtain the same results. Only the Pauli exclusion principle is at work. 
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APPENDIX A: COMMUTATOR [N, ^ (x)ip(x)] 
We want to show that 



with 



N = J2 d 3 k[4(k)c r (k) + dl(k)d r (k)\ . 

r J 

We use the following relation (F stands for fermion): 



Let's recall the expressions of the spinor fields: 



By using the anti-commutation relations 

{c s (k),4(p}} = 5 sr 5 3 (k-p) {d s (k),dt(p)} = S sr S 3 (k-p) 

and all other anti-commutators vanishing, we find that 

{^(f),c r (£)} = 

{^ a (f),4(fc)} = 

{ip a (x),d r (k)} = 
{^ a (x),4(k)}- 



1 Pm 



i m 
1 / m 
1 frn 



(Al) 
(A2) 



[FA F 3 F^ =F 1 [F 2 , F 3 F 4 ] + [F u F 3 F 4 ]F 2 (A3) 
=F 1 {F 2 , F 3 }F 4 - FiF 3 {F 2 , F 4 } + {F 1; F 3 }F 4 F 2 - F 3 {F 1 , F 4 }F 2 . (A4) 



(A5) 
(A6) 



(A7) 



(2vr) 3 V^ 



ul{k)e- i% - s 


{ip a (x),c r (k)} = 





(A8) 


U a r(k)e 


mx),ci(k)} = 





(A9) 


■Var(k)e-^ £ 


{4>t(x),d r (k)} = 





(A10) 




{ij a (x),dt(k)} = 





(All) 
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so that 



[V>i(£)Va(£),E / d 3 kct(k) Cr (k)] = (A12) 
r J 

E / rf 3 ^W(5){^(5),4(^)}cr(^) -cJ(fc)M(f) > c P (fc)R(5)) = (A13) 

77^ E / ^f^^tf)^ 

W ,,r ^ V E P E k 

(A14) 

m 2 >^ /" d 3 pd 3 k t 



—- ^ E / -^=[ut(k)uMe^~ ky "ci(k)cM + 4(fc)^(p)e-^)-ct(fc)4(p)] . 

(A15) 

Since r, s, p and are dummy variables, we find that 

[V4(£)Va(£),E / ^4(^r(£)] = (A16) 

E / ■^S=[«!(p)«r(*)e i ^- 2 d.(p)c r (fc)]- (A17) 



m 2 



(2tt; 



m 



3 _ 



JlE / ■^^WC^ti.^e^VcfC^c^] . (A18) 
In the same way, we obtain 

[^MMx), E / ^4$*®] = (A19) 

E / rf*fc(-^(5)d!:(*){^a(5), + 4(*)K(*W«(*)) = - (A20) 



(A21) 



(A22) 
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This can be simplified to 

[i;l(x)Mx),J2 I d 3 kdt(k)d r (k)] = (A23) 

r J 

jfv, E / 4=H(p)^(^)^^ ) ^(p)4(fc)]+ (A24) 
fS^E / ^^[vmu.W^^drfocM . (A25) 

W s ,r J v e p e u 

Putting the two results together, we get 

[^(£)VU£),iV]= (A26) 
|^E / ^=[«!(P H )«r(*)e- <(M) ^:(pl4(*)]+ (A27) 

I^E / 4ftK t (^)^(p1e^^)-X(^)c s (p1] , (A28) 

which is not equal to zero, even if we think about fields as distributions. If we start from 
the state <4(po)4(Po)|0), it is clear that there are well-behaved functions / such that 

(0| j d 3 xf(x)[^(x)^(x), N]\dl(p o )cl(p o )\0) + . (A29) 
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